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ABSTRACT 

We present results of fully self-consistent iV-body simulations of the motion of 
four globular clusters moving in the inner region of their parent galaxy. With 
regard to previous simplified simulations, we confirm merging and formation of 
an almost steady nuclear cluster, in a slightly shorter time. The projected surface 
density profile shows strong similarity to that of resolved galactic nuclei. This 
similarity reflects also in the velocity dispersion profile which exhibits a central 
colder component as observed in many nucleated galaxies. 

Key words: stellar dynamics - methods: numerical - galaxies: kinematics and 
dynamics - globular clusters: general 



1 INTRODUCTION 

The motion of massive globular clusters (GCs) in 
a galaxy is determined by both the smooth general 
potential and the graininess of the stellar field. This 
latter, fluctuating, component acts as a decelerating 
mechanism: the so called dynamical friction (hereafter 
DFR). The role of DFR, which dissipates cluster orbital 
energy and angular momentum, has been found to be 
crucial in determining the time evolution of globular 
cluster orbits in a galaxy, especially for orbits plung- 
ing in the inner, high density galactic regions (e.g. 
ICapuzzo-Dolcetta Sz Vicaril 1200a ). The decay time may, 
indeed, be short enough for massive clusters to limit 
their motion to the inner galactic regions. However, 
any general astrophysical consideration regarding dy- 
namical friction cannot be based just on approximated 
evaluations based on Chandrasekhar-like formulas, even 
in their gener alizations apt to treat axisymmetric or 
triaxial cases JPesce. Capuzzo-Dolcetta fc Vietril Il992 
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Several effects are indeed neglected by the analytical 
approach. For example, it is not considered the change 
in magnitude of the DFR due to the increase of the spa- 
tial extension of the cluster caused by the formation, and 
subsequent 'expansion', of its tidal tails, and, moreover, 
how the DFR acts on the different parts of the cluster 
(e.g., it could be stronger on its core than on the tails). 
Recent numerical experiments have tried to shed light 
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on this aspect finding that stars stripped from the clus- 
ter by the field, but still close enough to the system, 
continu e to contribute to the mass of the decelerating 
system l|Fellhauer fc Linll2007n . and that, in general, the 
real DFR effect can be stronger than that estimated by 
the usual Chandrasekhar formula (probably because of 
the further friction due to tidal effects, see Fuj ii et al.l 
120071 '). Another unexplored question is how the gravita- 
tional feedback on the very inner part of the galaxy in- 
fluences the DFR strength; this may be important when, 
during the final stages of orbital decay, the GC orbit gets 
very near to the centre so to enclose a galactic mass com- 
parable with that of the cluster itself. 

Clarifying the role of the above-mentioned dynamical 
effects is important also in the attempt to understand the 
mechanisms leading to the formation of Nuclear Clusters 
(NCs). In fact, it can help in discriminating the various 
scenarios proposed, especially in supporting the validity 
of multiple merging of s maller sub-systems as a forma- 
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tion m echanism (se e , e.g. 
| 2002| : iBekki et al.l 12004 
120081 ). At this regard, we cite the recent observational 
evidence of the existence of a very young and massive 
star cluster in NGC 2139, located 2" offset from the 
kinematical centre of the host galaxy, suggesting a for- 
mation that was independent and non coeval with that 
of the galactic bulge and, moreover, the observed envi- 
ronment is such that the system can decay to the centre 
in a time so short to keep intact its structure and be- 
come w hat is normally ca l led NC (Anderse n et al 1 120081 . 
see also iBoker et al]|2002l : IWalcher et alj|2006j h 

Reliable indications would follow by straightforward, 
direct A^-body integrations, which can be enormously 
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time consuming when treating self-consistently the mo- 
tion of a GC in a dense galactic environment. In this 
Letter we present the preliminary results of a fully self- 
consistent iV-body simulation concerning the close inter- 
action of a sample of four massive GCs in the central re- 
gion of a galaxy. Both the clusters and the galaxy are rep- 
resented by mutually interacting particles, thus including 
self-consistently both DFR and tidal interactions. This 
study can give depeer insight into the decay and merg- 
ing of GCs in galactic nuclear regions, a scenario first 
tackled by semi-analitical ap proaches (jTremaine et al.l 
Il975l ; ICapuzzo-Dolcettalll993h and then pursued by one 
full ./V-body experiment though w ith a resolution much 
lower than that presented here ()Oh fc Linl l2000l . who 
used N = 500 'particles' for each GC and N ~ 10 4 
to represent the galaxy). Simulations presented here 
are useful to test and validate pr evious recent results 
l|Capuzzo-Dolcetta fc Miocchill20Qg| , hereafter CM08) ob- 
tained in a simplified scheme in which the clusters are N- 
body objects moving in a fixed analytical galactic poten- 
tial and subjected to a dynam i cal fri ction braking given 
analitically by the IPesce et al . (1992) formula. 



2 SIMULATIONS 

Our iV-body simulation has been carried out em- 
ploying a parallel tree-code with a leap-frog time 
integrator using individual and variable time steps 
l|Miocchi fc Capuzzo-Dolcettal l2002h . Th e simu lation is 
set in a way similar to that described in CMOS, but for 
the role of the external galactic field which was, there, 
represented as a triaxial, analytical potential while, here, 
it is self-consistently (iV-body) represented starting from 
initial conditions sampling a given equilibrium density 
profile. We resume briefly the main characteristics of our 
modelization. We choose to follow four massive GCs as 
templates of a multi-merger among stellar systems de- 
cayed, conserving their individuality, in the inner galaxy 
region. A quick orbital decay induced by DFR is due 
to initial large values for the total mass of the clusters, 
whose resistance to external tidal disturbance was guar- 
anteed by their sufficient initial compactness. The GC 
sample corresp onds to the four densest clusters dealt 
with in ICMOSl and, moreover, every GC is composed by 
N — 2.5 x 10 5 equal mass stars. See Tab.[T]for the initial 
values of their structural parameters. 

The stellar bulge where GCs move is here repre- 
sented as another iV-body system, with Nb = 512, 000 
particles, sampling a Plummer model such to give the 
bodies initial positions and velocities according to a Mon- 
tecarlo representation of its equilibrium distribution func- 
tion (oc [— E] 7 ^ 2 with E being the particle total energy 
per unit mass). The Plummer density profile 



Pb(r) = p h0 [l + (r/r hc y 



-5/2 



(1) 



where pbo = Sy/SMb /47rr|j c is the central density, has two 
scale parameters, the core mass (Mb) and the core radius 
( r bc) which are used, in the following, as mass and length 
units, respectively. Consequently, time, velocity and den- 
sity will be measured in units of the galactic crossing 



time t b = {r'iJGMb) l/2 , of v h = (GM b /r bc ) 1/2 and of 
the central density pbo, respectively. Due to the infinite 
extension of the Plummer sphere, it has to be truncated 
at a distance r cut large enough to guarantee the stabil- 
ity of the sampling iV-body representation. This is done 
by the choice of r cut = 12.2rb c so as to contain 99% of 
the total mass. In order to have nearly t he same galac- 
tic mean gravitational field of that in the CMOS model, 
we set the Mb and rbc values as equal to those of the 
corresponding parame ters of the galactic density distri- 
bution used in ICMQq . even though this gives a slightly 
larger pbo- The four clusters [(a), (b), (c) and (d)] start 
th eir evo lution from the same initial conditions chosen 
111 CMOS. As clusters reference centres, we chose their 
centre -of-densities (CDs) as defined in lCasertano fc Hull 
( 1985). Finally, note that the simulation results are scale- 
invariant: any quantity can be re-scaled by fixing arbi- 
trarily the values of Mb and rbc in physical units. 



3 RESULTS 

The rapidity of the merger is evident in Fig. [T] as well 
as in the upper panel of Fig. [3] the merger process is 
completed at about t m ~ 17 tb, when the Lagrangian radii 
of the four GCs, seen as a whole system, flatten to a quasi- 
stable configuration (see Fig.[2j upper panel). In physical 
units: 
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100 pc 
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-1/2 



Myr. 



(2) 



Of course this 'merging' time depends much on the or- 
bital initial conditions of the progenitor clusters and mea- 
sures the time of the merger since the time in which 
the GCs are already confined within the galactic core 
region. Nevertheless, given that DFR has been convinc- 
ingly shown to be an efficient mechanism to drag GCs 
in the very inner regions of a ga laxy within a time much 
shorter than a Hubble time (see | Capuzzo-Dolcetta|[l993l . 
Capuzzo- Dolcetta fc Vicaril [2003 and the discussion in 
Miocchi et al.ll2006l ). this result supports the hypothesis 
stating that a GCs 'sedimentation' can take place at the 
galactic centre well within the galaxy life-time. 

It is interesting noting (see the upper panel of Fig. [3} , 
that the CD of the merger remnant oscillates within 
0.02rb c from the galactic CD. If, e.g., rbc ~ 200 pc this 
displacement corresponds to ~ 4 pc, i.e. a separation that 
would be hard to be resolved (corresponding to ~ 0.05" 
at the Virgo cluster distance), so the NC would appear 
substantially centered. The off-centered NC position is 
due to the gravitational feedback between the very cen- 
tral part of the galaxy and the clusters, which naturally 
derives from the self-consistent nature of the model. Such 
an interaction reflects into the 'perturbed' behaviour of 
the 1% Lagrangian radius (below ~ 0.2rb c ) of the galaxy 
that is evident from Fig. [2] (lower panel) . 

The lower panel of Fig. [3] indicates that a first merg- 
ing event between 2 clusters [(b) and (d)] occurs rather 
early (at t < 5), when their centres are still ~ 0.2rbc 
far from the galactic centre. The fact that such a bi- 
nary coalescence takes place before the complete orbital 
decay, suggests that these events can occur even offset 
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Figure 1. Time sequence of the four GCs merger. Field stars are not shown. 



Table 1. Parameters list for the initial cluster models, expressed in galactic units. Reported 
are: the GC mass (M), the limiting radius (rt), the King radius (r c ), the King concentration 
parameter (c), the half-mass radius (r^), the central density (po)> the half-mass crossing time, 
t c h = [rJ/fGM)] 1 / 2 , and the King velocity parameter (<r K ). 
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Figure 2. Evolution of the Lagrangian radii of the four GCs as 
a whole (upper panel) and of the galaxy (lower panel) averaged 
over a time window = 3<b- They correspond to the radii of 
the sphere centered at the system CD and enclosing a given 
fraction of the total mass (as indicated). 



Figure 3. Upper panel: time evolution of the distance of the 
cluster CD from the galactic centre (rcr>). Lower panel: time 
behaviour of the distance between each cluster CD and its 
first-neighbour. All the behaviours have been averaged over a 
time window = 3tb- Solid lines refer to cluster (a); dotted: (b); 
short-dashed: (c); long-dashed: (d). 



from the galactic centre, in spite of the strong tidal field. 
This 'pre-decay' coalescence was also found in the low- 
resolution self-consistent simulation bv lOh fc Linl (2000) 
(their model lb), though in their case this phenomenon 
takes place much farther from the galactic centre. 

By comparing the merger time with that of the 
|CM08| simulation, it is found that here the orbit al decay 
is ~ 1.6 times faste r than that given by the use (in[CM08) 
of the generalized l|Pesce e t al. 1992) Chandrasekhar for- 
mula 

This means that the DFR efficiency is greatly en- 
hanced in the self-consistent model. Indeed, in Fig. [4] we 



see that, at least for t < 20, the clusters CDs are always 
clo ser to t he centre of the galaxy respect to the case of 
the ICM08I simulation. For later times the comparison is 
no longer reliable, because of the very small values of rcD ■ 
In Fig. [S] the total surface density profile is plotted 
for the last configuration and the typical appearance of a 
nucleated galaxy central profile comes clearly out. As an 
example, the s imilarity with the VCC 1871 profile in the 
Virgo cluster |C6te et all 120061 ) is particularly evident. 
The central 'over- density' associated to the NC presence 
is, howev er, less pronounced with respect to that obtained 
m ICMOSl and this is due to the lower concentration of the 
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Figure 4. Time evolution of the ratio between the distance of 
the cluster CD fr om the galactic centre and the same quantity 
evaluated in the CM08j simulation (rcD,o)i averaged over a 
time window = 3tb • Line symbols are the same as in Fig. [3] 
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Figure 5. Projected surface density profile of the whole sys- 
tem (galaxy plus NC) in the final configuration (solid line) and, 
for the sake of comparison, of the galaxy only (long dashed 
line). 



NC in comparison with the last configuration in ICMOSl 
case. Perhaps, the lower NC central density may be a con- 
sequence of the too large collisional effect generated by 
the interaction of the clusters with the 'particles' repre- 
senting the very inner region of the galaxy that, because 
of the limited resolution, are probably too massive. Only 
higher resolution simulations of the very last stages of the 
orbital decay and merging would allow to establish the 
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Figure 6. Projected radial behaviour of the velocity disper- 
sion in the last system configuration. Line symbols are as in 
Fig.E 



actual role of this effect (Capuzzo-Dolcetta & Miocchi, in 
preparation) . 

An important result concerns the peculiar behaviour 
of the gl obal vel ocity dispersion profile. Confirming the 
finding in CMOS, the velocity dispersion of the whole sys- 
tem (galaxy plus NC) shows a clear decrea se towards the 
centre (see Fig. |6}. As discussed in lCMOct this is a clear 
sign of the presence of two kinematically distinct sys- 
tems that are relaxed into a different dynamical status. 
Such a peculiar feature has been actually found in most 
of the Virgo cluster nucleated dwarf ellipticals observed 
bv lGeha. Guhathakurta fc van der M arel ( 20021. see their 
Fig. 5). The decreasing trend to the centre of the veloc- 
ity dispersion is also consistent with the solution of the 
Jeans e quations for a sample of NCs observed in late- type 
spirals dWalcher et al.ll2005l ). Finally, either a flattening 
or a slight centr al decrease are found and discussed in 
lOh fc Lin! l|2000l ). too. 



4 SUMMARY 

Massive globular clusters suffer of dynamical friction 
braking such to be limited to move in the inner part of 
the host galaxy. The interaction of decayed globular clus- 
ters among themselves and with the external tidal field 
may lead to subsequent merger events, whose modes are 
not easily understood with out det ailed ./V-body simula- 
tions. In a previous paper |CM08) we followed numeri- 
cally the interaction and merging of four globular clusters 
each represented as an TV = 250, 000 stars system mov- 
ing in an analitical triaxial potential , taking dynamica l 
friction into account by means of the lPesce etafl (|l992h 

formula. 

The present paper has the aim to check the ICMOSl 
results via a full self-consistent JV-body modelization, i.e. 
sampling also the galactic environment by a large num- 
ber (= 512, 000) of interacting particles. All the main 
features of the me rging e vent and of the 'super cluster' 
remnant found by ICMOSl are qualitatively confirmed by 
this, more detailed, simulation, but for a stronger orbital 
decay (two times faster) due to collective effects and for 
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a less centrally concentrated merger product. This latter 
result, however, needs to be checked by higher resolution 
simulations which are also necessary to study reliably the 
long-term evolution and stability of the formed NC. 
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